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ABSTRACT 

We determine the distance to the open cluster NGC 2264 using a statistical analysis of cluster 
member inclinations. We derive distance-dependent values of s'mi (where i is the inclination angle) 
for 97 stars in NGC 2264 from the rotation periods, luminosities, effective temperatures, and projected 
equatorial rotation velocities, v sin i, measured for these stars. We have measured 96 of the v sin i values 
in our sample by analyzing high-resolution spectra with a cross-correlation technique. We model 
the observed distribution of sini for the cluster by assuming that member stars have random axial 
orientations and by adopting prescriptions for the measurement errors in our sample. By adjusting 
the distance assumed in the observed sini distribution until it matches the modeled distribution, 
we obtain a best-fit distance for the cluster. We find the data to be consistent with a distance to 
NGC 2264 of 913 pc. Quantitative tests of our analysis reveals uncertainties of 40 and 110 pc due to 
sampling and systematic effects, respectively. This distance estimate suggests a revised age for the 
cluster of ~1.5 Myrs, although more detailed investigations of the full cluster membership are required 
to draw strong conclusions. 



Subject headings: methods: statistical — stars: formation 
open clusters: NGC 2264 



stars: rotation — stars: distances 



1. INTRODUCTION 

NGC 2264 is an open cluster in the Monoceros OBI 
association containing a large population of young stars. 
The cluster has been the focus of a number of st udies of 
early stellar evolution dating back to the work of lHerbigl 
(1954). NGC 2264 is an ideal target for such studies 
because of its low line-o f-sight extinction and minimal 
optical nebular emission (jPark et al.l l2000). In addition, 
most background stars are obscured from the cl uster by 
the p resence of a large molecular cloud complex (Hcrbig 

QUI). 

Previous distance determinations for NGC 2264 have 
favored a distance of 800 pc, although there h as been 
consid erable spread around this value: 800 pc by|Walkcr 
(fl956h. 760 pc b v [Su]i7et~all (EMZl) , 760-950 pc by 
iFlaccomio et al.1 <u999ft . 750 pc by iMavne fc Navlorl 
(|2008f) : see lDahml (|2008ft for a complete summary of pre- 
vious distance measurements. These distance estimates 
have been primarily derived by fitting empirically or the- 
oretically defined main sequences to the location of high 
mass (B and A type) stars in the HR diagram. Although 
extensive Stromg ren narrowband photometry is availa ble 
in the literature (|Strom et al.lll971tlPerez et i al.lll989f) . it 
has not yet bee n used to refine distanc es to these early 
type stars (e.g., IAnthonv-Twarodll982fl . 

An improved distance estimate to NGC 2264 would 
lower the uncertainties in the luminosities derived for 
cluster members, which in turn constrain the cluster's 
age. As one of the premiere laboratories for studying 
star formation in the Milky Way, a better estimate of 
the age of NGC 2264 would further our understanding of 
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processes relevant to early stellar evolution such as angu- 
lar momentum transfer and the lifetime of circumstellar 
disks. 

In this paper, we determine the distance to NGC 2264 
using a statistical technique that relics on measured pro- 
jected rotation velocities, rotation periods, luminosities 
and effective temperatures of the low mass K & M type 
cluster memb e rs. T he technique was first developed by 
I Hendry et al.l (|1993l ) and has subsequently been used to 
find distances to the Pleiades, the Tau rus star forming re- 
gion, and the Orion Nebula Cluster (jO'Dell et al.l ll994l: 
iPreibisch fc SnTrtllfll9^lJeffriesl[2W^ . This method has 
the advantage of being nearly independent of stellar evo- 
lutionary models. 



1.1. The Method 

In brief, we first measure the projected rotational ve- 
locities of cluster members, vsmi (where v is the tan- 
gential velocity of the stellar surface at the equator and 
i is the inclination of the stellar rotational axis on the 
sky such that i = 90° implies an edge on orientation and 
i = 0° implies a pole on orientation) from existing high- 
resolu tion spectra of NGC 2264 members (Fure sz et al.l 
2006). An effective temperature, T e //, is estimated for 
each star from either its spectral type or dereddened 
photometry. Luminosities, L, are estimated for cluster 
members from measured magnitudes by assuming a nom- 
inal value for the cluster distance, a prescription for the 
cluster reddening, and a standard bolometric correction. 
Stellar radii are then calculated from the estimated lu- 
minosities and effective temperatures using the Stefan- 
Boltzmann relation. The final data needed for the dis- 
tance determination are rotation periods obtained from 
fits to periodic variations in the stellar light curves. 

Bringing all of these data together, we calculate sini 
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for each star as 



sin! 



P- (v smi) 
2irR D 



(1) 



where P is the measured rotation period and Rd is the 
stellar radius (the subscript indicates that the measured 
radius is dependent on the stellar luminosity, which de- 
pends on the assumed cluster distance, D). 

Eq. Q] allows us to generate an observed distribution of 
sin i that is dependent on the adopted cluster distance. 
We then model the distribution of sini assuming that 
the rotational axes of stars in the cluster are randomly 
oriented. By scaling the input cluster distance so that 
the observed sin i distribution matches the predicted dis- 
tribution, we obtain an estimate for the true cluster dis- 
tance. 

We begin in <|2] with a description of the data used 
to determine the distance to NGC 2264; ^describes our 
measurement of projected rotation velocities by applying 
a cross-correlation routine to high resolution spectra of 
cluster members; ^describes the distance determination 
technique in detail; our results are discussed in i j4.61 

2. DATA 

As described above, our distance measurement relies 
on sini values for stars in NGC 2264. Four types of data 
are needed to calculate sini for an individual star using 
Eq. [T] the star's period, luminosity, effective tempera- 
ture and projected equatorial rotational velocity, usini. 
Period (i j2.1|) . effective temperature ( §2.2j) and luminos- 
ity (J^2l_data wereobtained from a catalog compiled 
by iRebull et alj (|2006l ) and are briefly described below. 
Most of the v sin i data, on the other hand, were calcu- 
lated from spectra ( §2.4|) using a cross-correlation tech- 
nique further discussed in 

2.1. Periods 

Rotation periods for pre-main sequence stars are de- 
termined by measuring periodic variations in the objects' 
brightness. These variations arise from the presence of 
large (~ 40° angular radius) starspots o n the surface s 
of these young, magnetically active stars (|Herbstl fl989f ). 
The transit of the starspot(s) as the star rotates dimin- 
ishes the observed stellar luminosity on the order of a 
few tenths of a magnitude. This variation is quite stable 
a nd can be use d to de rive precise periods. 

IRebull et al.1 (|2006h compiled periods meas ur ed fo r 
members of NGC 2264 bv IRebull et al.l (2002), 



iMakidon et~aT1 (|2004h and lLamm et all (|2005T ). The pe 
riod measurements range from roughly 0.5 to 29 days and 
are predicted to be accurate to roughly 5P/P ~ 1%. 

2.2. Effective Temperatures 

The catalog assembled by Re bull et al.1 (|2006l ) includes 
T e f f estimates derived from low-resolution spectral types 
and from the stars' optical colors using the V — I vs. 
T e ff relation presented by Hil lenbrand ()1997t ). Each 
star's V — I color was dereddened p rior t o this calcu- 
lation as described bv IRebull et alj (|2002f ): stars with 
spectra were dereddened so that their observed R—I col- 
ors matched the intrinsic colors of that spectral type on 
the zer o age ma i n seq uence (ZAMS) defin e d by iBesselH 
(|1991h . iLeggettl (|1992ft and iLeggett et all (|1998l ). The 



V— I colors for stars without spectra were dereddened as- 
suming the modal reddening of members with measured 
spectral types. While the use of an overall reddening for 
the cluster is less than ideal, out of a total of 97 stars for 
which we derive sini values, only 14 lack spectral types. 
Thus, larger errors associated with photometric based 
T e f f measurements likely have a negligible effect on our 
results. 

T Tauri stars almost certainly do not have a single 
photospheric temperature, however ; obse rvations of the 
weak T Tauri star V410 by iHerbstl (|1989D indicated the 
presence of two large, polar star spots with characteristic 
temperatures of 3100 K, in constrast to the star's 4400 
K. photosphere. Observations of larger ensembles of T 
Tauri stars indicate similar photospheric-spot tempera- 
ture di fferentials, and typical spot covering fractions of 
~10% dBouvier fc Bertout IH989 : 1 Johns-Krull & Gaffordl 
12001 . 

Somewhat counterintuitively, however, star-spots are 
unlikely to introduce large errors in a star's derived 
photospheric temperature, particularly for large spot- 
photosphere temperature differentials. The star's in- 
tegrated emission is dominated by non-spotted photo- 
spheric flux, as the cooler spots emit significantly less flux 
per unit area. O bservational confirm ation of this effect 
was provided by iFrasca et al.l (|2005l ). who constructed 
detailed star-spot models of RS CVn systems to repro- 
duce broadband photometric light curves and tempera- 
ture sensitive line ratios from temporally resolved, high 
resolution spectra of RS CVn system s. Whi le the best 
fit models identified by IFrasca et alj (|2005f ) possessed 
photosphere-spot temperature differentials of ~1000 K, 
the temperature sensitive spectroscopic line-ratios only 
departed from the photospheric value by ^150 K over 
the course of the obser vations. 

In their Appendix A, IFrasca et al.l (|2005f ) outline a for- 
malism for calculating the mean temperature observed 
from a star with a given spot covering fraction and 
photosphere-spot temperature differential. Using this 
formalism, we calculated the difference between a T 
Tauri star's true (non-spotted) photospheric tempera- 
ture and the temperature that would be measured from 
its integrated (spot + photsphere) spectrum. For typi- 
cal T Tauri star parameters (T p/lot =4000 K; T spot = 3200 
K; f sp ot = 0.1), the temperature measured from the com- 
bined spot+photosphere I-band spectrum would be 3992 
K, only 8 K different from the 'true' photospheric tem- 
perature. 

Given the small size of this effect, we assume the un- 
certainties in our derived temperatures are dominated by 
errors in the observed colors, spectral types, and redden- 
ing corrections. We expect the temperature uncertainties 
introduced by these effects are roughly 5T e j f/T e j / « 5%. 

2.3. Luminosities 

IRebull et al.l (|2006l ) calculated luminosities for cluster 
members from dereddened photometry, assuming a dis- 
tance to the cluster of 760 pc. Stars without previously 
determined spectral types were dereddened by the modal 
reddening value determined for stars with spectral types. 
Not taking into account the uncertainty associated with 
the distance assumption, the observed luminosities are 
likely accurate to within SL/L 35%. This estimate in- 
cludes contributions from uncertainties in spectral type, 
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exti nction values, an d uncertainty due to source variabil- 
ity (|Hartmannll2001h . 

We note that the luminosity we calculate for each 
star from its measured photospheric temperature will 
likely be a slight over-estimate, as some fraction of the 
star's surface will be covered with cooler, less lumi- 
no us star-spots . Aga in using the formalism presented 
by iFrasca et alj (|2005l ) in their Appendix A, we calcu- 
lated the expected scale of this effect. For our typi- 
cal T Tauri star, the bolometric luminosity calculated 
assuming a single measured photospheric temperature 
(ie, T Q fc s =3992 K) is 5% larger than the true luminosity 
assuming the correct spot-photosphere differential and 
covering fraction. As the scale of this effect is well within 
our errors, and as we lack a robust characterization of the 
spot properties for the stars in our sample, we simply 
adopt the luminosity implied by the single temperature 
associated with each star's spectral type. 

2.4. Echelle Spectra 

Spectra for 923 stars in and aroun d NGC 2264 were ob - 
tained from the sa mple observed by iFuresz 
using Hectochelle (jSzentgvorgyil l2006h . a multiobject 
echelle spectrograph located on the 6.5 m Multiple Mir- 
ror Telescope (MMT). The spectra cover a wavelength 
range of 6450 A to 6650 A and have a resolution of 
R ~ 34, 000, providing velocity resolution of ~ 9 km/s. 
A typical spectrum obtained for a star in NGC 2264 is 
shown in Fig. [TJ 

3. PROJECTED ROTATIONAL VELOCITIES 
3.1. The Correlation Technique 

Because few stars in our initial catalog had previously 
measured Dsini's, it was necessa r y to e xtract this data 
from the spectra of IFuresz et alj (|2006l ). To make this 
measurement, we used a cr oss-correlation t e chniq ue sim- 
ilar to that developed by iTonrv fc Davisl (jl979l ). The 
correlation parameter, C(vr), between a target (i.e. a 
cluster member) and a well-matched template spectrum 
is determined as a function of the radial velocity of 
the template spectrum, vr. C(vr) is obtained by in- 
verse Fourier transforming the product of the discrete 
Fourier transforms of the target and template spectra 
(|Hartmann et al.l fl9861 . To generate C(vr), we have de- 
veloped a custom cross-correlation routine in IDL based 
on a heavil y modified version of a rout ine originally de- 
veloped bv 1 White fc Hillenbrand! (|2004D . 

The location of the peak of C(vr) corresponds to the 
value of vr for which the target and template spectra are 
best matched. Thus, the correlation procedure provides 
us with a measure of vr for stars in our sample. Further- 
more, because rotational broadening of the target star's 
spectral lines results in a broader C(vr) peak, the width 
of the C{vr) peak, ac, provides a measure of wsini for 
the target star. The correlation function C(vr) resulting 
from cross-correlating one of our target spectra with a 
template spectrum is shown in Fig. [2j 

The accuracy of the kinematic properties measured by 
the cross-correlation technique depends on the quality 
of the agreement between the strength and shape of the 
spectral features present in the template and target spec- 
tra. To ensure that the target and template spectra are 
reasonably well matched, we correlate each target spec- 
trum against a grid of template spectra covering a range 



of temperatures and surface gravities (see ^3.2(1 . The 
best matching template is then used to derive v sin i for 
that target. The degree of 'matching' between a template 
spectrum and a target spectrum is quantified with a sta- 
tistical quantity, R, defined as the ratio of the height of 
the maximum peak in C(vr) to the root -mean-square of 
the a ntisymmetric component of C(vr) (|Tonrv fc Davisl 
119791 ). In an idealized scenario in which the target spec- 
trum has no noise, the correlation function would be per- 
fectly symmetric around some particular value of vr. R 
therefore measures the strength of the correlation func- 
tion peak against the noise in C(vr). Thus, the tem- 
plate that produces the largest R value when correlated 
against a particular target can be said to be the optimal 
template for that target. 

3.2. The Template Grid 

We cross-correlat ed our targets again st a set of tem- 
plates calculated bv lCoelho et al.1 (|2005f ). The template 
grid covers temperatures ranging from 3000 K to 7000 K 
in increments of 250 K and log g (surface gravity) values 
ranging from 0.0 to 5.0 in increments of 0.5. Because the 
strength of the Ha line at 6562 A exhibits strong depen- 
dence on accretion that could potentially interfere with 
our ability to derive wsini values, each spectrum was di- 
vided into two regions bracketing Ha (a 'blue' region at 
wavelengths below 6545 A and a 'red' region at wave- 
lengths greater than 6585 A). Both regions were then 
correlated separately with the template spectra. Subse- 
quent analysis revealed that our ability to derive accu- 
rate v sin i values was greater for the cross-correlations 
performed on the red region. This effect was likely due 
to the fact that the C(vr) peaks tended to have a more 
Gaussian shape for the red cross-correlations. Thus, in 
future discussions R refers to the value of R calculated 
for the red region. 

The relationship between ac and the target v sin i 
was calibrated separately for each template spectrum by 
cross-correlating the template spectra with synthetic tar- 
get spectra for which vsini is known. The synthetic 
target spectra were created by introducing artificial ro- 
tational broadening into the template spectrum and de- 
grading the spectral resolution and sampling to match 
the Hectochelle data. We performed this test for each 
template with synthetic targets whose usini's ranged 
from to 99 km/s in steps of 3 km/s. The best fit- 
ting template for the example target spectrum in Fig. 
[1] is also shown in that figure, modified to reflect the 
corresponding vr and v sin i as determined through our 
cross-correlation procedure. 

3.3. Parameters Derived from Cross-Correlation 

Once the best matching template spectrum for a tar- 
get has been found, the target's ushn is measured by 
using the width of the cross-correlation peak as an in- 
put for a fifth-order polynomial fit to the ac-vsini re- 
lation previously determined for that template. Valid 
vs'mi values can only be derived when the ac for the 
target is within the range of ac calculated for the tem- 
plate. In total, valid usini values were derived for 489 
stars. Fig. [3] shows the distribution of usku values for 
stars in our sample. The absence of stars with vsini < 
6 km/s likely represents a combination of a true lack of 
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Fig. 1. — A spectrum of a typical star in our sample (solid line). Also plotted (dashed line) is the best fitting templa te sp ectrum for this 
star (see §3j2}. The template spectrum has been shifted to lower flux values for visual clarity. As noted in the text ( i|3,3|l the Ha region 
(from 6545 A to 6585 A) was excluded from the cross-correlation procedure because of its strong variability between spectra. 
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Fig. 2. — A cross-correlation function resulting from the corre- 
lation of a target spectrum with the corresponding best fit tem- 
plate model. The cross-correlation function shows a clear peak at 
roughly 10 km/s; the width of the function formally implies a rota- 
tion velocity of about 3.4 km/s, below our threshold for confident 
vs'mi detections. A Gaussian fit to the cross-correlation function 
is shown as a dashed line. The above example represents a better 
than average cross-correlation result, with an R value of approxi- 
mately 44. 

stars with low v sin i, as well as the resolution limit of the 
Hectochelle spectrograph (9 km/s). The finite spectral 
resolution and standard measurement uncertainties pro- 
duce cross-correlation widths for low v sin i stars that are 
narrower than the minimum ac values in the template 
calibrations. Our inability to separate the instrumental 
bias from the true statistics at low v sin i motivated our 
decision to impose a cutoff in v sin i (see £|4.2p . 

3.4. Uncertainties in usini 

Rigorous uncertainties were calculated for the v sin i 
values derived from the cross-correlation procedure by 
running a Monte Carlo test on a set of artificial spec- 
tra. The synthetic spectra were rotationally broadened 
to arbitrary v sin i and then matched to the resolution 
and sampling of the Hectochelle data. Gaussian noise 
was also added to the synthetic targets to simulate the 
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Fig. 3. — The distribution of derived values of vs'mi. As can be 
seen from the plot, the observed frequency of stars with a particular 
vs'mi declines at vsini < 6 km/s as a result of an instrumental 
detection threshold as well as a true lack of stars with very low 
v sin i. 



effects of a finite signal to noise ratio in the actual data. 
The 1-cr amplitude of the noise was one fourth the mean 
signal level, characteristic of the lowest quality data in 
the Hectochelle observations. Since the value of vsini is 
known for the synthetic spectra, the difference between 
the calculated v sin i and the true v sin i for the synthetic 
spectra provides a measure of the uncertainties inherent 
in our cross correlation procedure. 

Uncertainties in v sin i were determined as a function of 
R for each synthetic target. Fig. 2] shows the deviation, 
Avsmi, between the input and calculated value of vsini 
for each synthetic target as a function of R. As should 
be expected, Av sin i decreases with increasing R since 
higher R corresponds to a be tter match between the tar- 
get and template. Following lHartmann et ail (|1986f ). we 
quantify the Av sin i vs. R relation by fitting a curve of 
the form a/(b + R) to the la width of Aw sin i. We find 
that the best fit is Av sin % = (19.7 km/s)/(R-2.74). We 
apply this relation to the R values measured from actual 
target spectra to quantize our errors in measured values 
of vsini. 
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Fig. 4. — Error in the vsini values calculated by our cross- 
correlation routine as a function of R. Each data point represents 
a synthetic target that has been rotationally broadened to some 
v sin i (true) and then correlated against the template grid to find 
the best matching template, v sin i (calculated) was then deter- 
mined from the relationship between vsini and the width of the 
correlation function peak. The solid line shows a fit to the 1-a 
width of the distribution, while the dashed line represents a similar 
fit to the 3-<r width of the distribution (see text for fit parameters) . 

4. DISTANCE DETERMINATION 

The combination of period, luminosity, effective tem- 
perature and v sin i data allows us to calculate sin i for 
each of the stars in our sample. Since sin i depends on the 
absolute luminosities of stars in our sample, the observed 
distribution of sin i values is inherently dependent on the 
assumed cluster distance. We use a Monte Carlo routine 
to produce model sini distributions that assume ran- 
domly oriented stellar rotation axes and that incorporate 
the effects of observational uncertainties. The distance 
to NGC 2264 is then constrained by comparing model 
sini distributions to observed distributions with differ- 
ent assumed cluster distances. The distance at which 
the modeled and observed distributions of sin i agree is 
our best fit distance. 

4.1. Calculating sini 

Assuming that the measured period, P, of each star is 
equal to its rotation period at the equator, we have 



sm i - 



P- (vsini) 



(2) 



2ttP ' 

where R is the stellar radius. P and vsini are direct 
observables, but R must be inferred from each star's lu- 
minosity, L, and effective temperature, T e ff, using the 
Stefan-Boltzmann relation. L must itself be estimated 
from each star's extinction-corrected bolometric magni- 
tude (uiboi) and an adopted distance to the cluster: 

_ ln -0.4(m bo! -51og_D+5) 



L = 10 



(3) 

Thus, we can write an expression for sini that is explic- 
itly dependent on distance: 



sini = 



P(v sin i)T e 2 



Since luminosities were previously calculated by 
iRebull et al.l (III) assuming a distance D — 760 pc, 
it is convenient to re- write Eq. 0] as 



P(vsmi)T^ 



10-0.2(m 6oi -51ogD+5) 



(4) 



(5) 
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where Lp is t he lum inosity that has been calculated by 
IRebull et alj (2006) for the stars in our sample. Eq. [5] 
allows us to calculate the sini distribution of stars in 
NGC 2264 from P, usini, T e // and Lq data, as well as 
an initial estimate of the cluster distance, Dq. 

4.2. Selection of Sample Stars 

Beginning with our sample of 923 stars for which we 
have Hectochelle spectra, we use a series of quality cuts 
to identify those stars that will allow us to produce an 
unbiased estimate of the distance to NGC 2264. We first 
restrict our analysis to the 489 stars for which we are 
able to measure v sin i using our cross-correlation rou- 
tine. We further restrict our sample to those stars with 
a high-likelihood of being bona-fide members of NGC 
2264. Since the cluster is a coherent kinematic system, 
imposing a radial ve locity cut ( 10 km/s < vr < 30 km/s) 
(|Furesz et al J 12006) on our spectroscopic sample identi- 
fies 273 likely members. Of these, 130 have the ancillary 
measurements (P, P e //, and Lr) necessary to compute 
sini. 

IRebull et al.l (|2006l ) identified a number of stars as 
likely members of NGC 2264 based on their positions in 
the sky and their location in color-magnitude space. We 
add one star to our sin i sample fro m the catalog o f NGC 
2264 cluster members compiled by IRebull et al.l (2006) ; 
this additional star is the only catalog member without 
a Hectochelle spectrum but with a previous v sin i mea- 
surement, as well as the other measurements necessary 
to estimate sini. Finally we restrict our sample to the 
97 stars with measured v sini values larger than 9 km/s, 
the minimum velocity resolution of our Hectochelle data. 
In the discussion that follows, we refer to this final sub- 
set of 97 stars as the 'distance sample'; for clarity we 
summarize the steps in its selection in Table [TJ 

Table [3] lists all the data used to calculate sin i for 
the 97 stars in the distance sample. The 'L& ; & T e ff 
source' column identifies the type of reddening correc- 
tion applied to each star's photometry in estimating its 
luminosity and temperature; stars with an 'S' have in- 
dividual reddening estimates based on observed spectral 
types, while stars labeled 'P' had photometry corrected 
assumin g the modal e xtinct ion derived for cluster mem- 
bers by IRebull et al.l ((2006). Finally, the star added to 
the dist ance sample using t he v sin i measurement cata- 
loged by IRebull et al.l (2006 ; lacks a corresponding v sin i 
error estimate: its entry in Table [3] lists 'None' in the 
vsini error column. 

Figure [5] shows the distribution of spectral types for 
stars in the distance sample. Most of the stars in the 
distance sample have spectral types in the range of K4 
to M3. 

4.3. Error Distributions for Observed Parameters 

Observational uncertainties affect the shape of the 
measured sini distribution. In order to obtain a reli- 
able distance estimate, we must account for these un- 
certainties in our modeled values of sini. To do this. 
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TABLE 1 

Summary of Distance Sample Selection 



Sample Subset 



Number of stars 



• Initial set of Furcsz et al. (2006) spectra 

• Reliable v sin i results 

• Radial velocity members w/ vsini 

• Ancillary data for sin i estimate, 
radial velocity member &c v sin i 

• v sin i > 9, ancillary data, 
& radial velocity member 



923 
489 
273 

130 

96 (+1 from Rebull et al. 2006) 
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Fig. 5. — The distribution of spectral types in the distance 
sample. The bulk of the stars in our sample have types between 
K4 and M3. 



we first estimate uncertainties in P, vsini, T e tf and L 
from the observational data; a Monte Carlo simulation 
is then used to incorporate these error distributions into 
the modeled value of sini. 

The rotation periods of T Tauri stars can be measured 
with high precision from their light curves. The errors 
associated with such measurements are usually on the 
order of 1%. In some cases, confounding factors such 
as aliasing in the light curves or the presence of multi- 
ple starspots can in crease these errors dramatically (e.g. 
iHerbst et al.l [2002V However, the number of cases in 
which these effects occur is typically small. Since the 
period errors are very small compared to the errors asso- 
ciated with other T Tauri measurements (i.e. luminosity 
and T e tf), we simply assume that fractional errors in pe- 
riod are normally distributed with a standard deviation 
of 1%, characteristic of the typical errors in pre-main 
sequence stellar period measurements. 

For the remaining variables it is possible to ascertain 
some measure of the actual errors from the observed 
data. Fractional v sin i errors for the distance sample 
were calculated using the relationship between R and 
S (vsini) given in £)3.4I The resulting error distribution 
is consistent with a normal distribution with <r=20%. 

Fractional errors in T e ff were determined by compar- 
ing values of T e ff calculated using individual spectral 
type-based reddening corrections or simply adopting the 
modal reddening for all cluster members. The T e // error 



distribution calculated using this prescrption is consis- 
tent with a normal distribution with a = 10%. We ex- 
pect, however, that calculating T e ff values by adopting 
the cluster's modal reddening is less accurate than deriv- 
ing reddenings from observed spectral types. Indeed, we 
note a clear relationship between the ST e ff value mea- 
sured for each star and its spectroscopic reddening esti- 
mate. This suggests that the difference between the two 
T e ff estimates is dominated by the errors introduced 
by adopting the modal reddening, and that the result- 
ing error distribution overestimates the actual errors as- 
sociated with the T e ff values derived using individual 
spectral type-based reddening corrections. 

Similarly, the errors in L are estimated from the dif- 
ference between luminosities calculated assuming a red- 
dening derived from each star's spectral type and those 
calculated assuming an overall reddening for the cluster. 
Wc find that the luminosity error distribution implied by 
these distinct L estimates lies within the bo unds of the 
norma l distribution with er=35% suggested by Hartmann 
(|2001| ) as characteristic of luminosity errors in pre-main 
sequence stars. 

4.4. Model sini Distributions 

Our technique for modeling the distribution of sini 
in NG C 2264 borrows heavily from iPreibisch fc Smith 
(|1997l ). We define the modeled value of sinz, (sini) m as 
follows: 



(sini) 



P (v smi) (T e f f )l^f<jjTrS 



Po (usin)o (T e// )g 



(6) 



Vl 

where the subscripted variables represent the actual val- 
ues independent of measurement uncertainties, and the 
non-subscripted variables represent the observed values 
including measurement uncertainties. Eq. [6] allows us to 
split the dependence of sini into two parts: the actual 
value of sin i (the term outside the brackets) and the con- 
tributions of measurement uncertainties (the term inside 
the brackets). 

The (sini) m distribution is generated from Eq. [5] us- 
ing a Monte Carlo routine. Each term of the form X/Xq 
(where X represents any of the variables P, usini, T e ff 
or L) is calculated by drawing randomly from the ap- 
propriate error distribution for our dataset, as described 
in £14.31 This process assumes that fractional errors in 
usini, T e ff and L are independent of the values of these 
variables. While this assumption is somewhat ques- 
tionable, the deviation from the true error distribution 
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is likely small. An example of a (sinz) m distribution 
that includes the effects of observational uncertainties is 
shown in Fig. [SJ 
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Fig. 6. — Modeled sini distribution assuming random axial ori- 
entations before (solid line) and after (dashed line) taking into 
account measurement uncertainties. The measurement uncertain- 
ties assumed in the generation of the model distribution are drawn 
randomly from the observed uncertainty distributions. We also 
restrict the model to tisini > 9 km/s (see i|4, 4.11 1. 



4.4.1. Distribution of vtrue 

In order to incorporate the effects of the v sin i cutoff 
adopted in £14.21 into our model it is necessary to assume 
some prescription for the distribution of the true equato- 
rial velocities, v true . Once such a prescription has been 
assumed, modeled wsini values can be calculated, allow- 
ing the model sample to be restricted in the same manner 
as the observational sample. 

We tested a variety of Vtrue distributions by combining 
randomly sampled v values with random axial orienta- 
tions and comparing the resultant set of modeled v sin i 
values with the observed v sin i values. The input Vtrue 
distribution was adjusted until the sampled population 
produced a satisfactory match with the observed v sin i 
values, as indicated by a two-sided Kolmogorov-Smirnov 
(KS) test. We find that modeling the v trU e distribution 
as an exponentially decaying function with a constant 
offset (P{v trU e) oc e~ a ' Vtrue + C) leads to a good match 
with the observed wsini distribution. Our best match- 
ing Vtrue model has a = 0.09 and C = 0.004. A KS test 
comparing the resultant modeled v sin i distribution with 
the observed distribution yields, on average, a probabil- 
ity of ~ 95% that our modeled v sin i distribution comes 
from the same underlying distribution as the observed 
wsini's. A flat distribution of Vtrue, on the other hand, 
can be rejected with a probability greater than 99.999%. 
See Fig. for a comparison of the assumed Vtrue distri- 
butions . This result is in agreement with that found by 
I Jeffries! (12001 . 

4.4.2. The Binary Correction 



Actual 

Exponential Model 
Flat 




40 60 
v sin(i) 



100 



Fig. 7. — Comparison of the usini distributions resulting from 
different assumed vtrue distributions. The exponential model as- 
sumes the probability of a star having a particular vtrue goes as 
an exponentially decaying function, plus a constant offset. The 
flat model assumes that all vtrue in the range considered (0 to 
100 km/s) are equally likely. As can be seen from the figure, the 
exponential model fits well to the data. 

Some unknown fraction, J5, of the stars included in our 
catalog are unresolved binary systems. For such systems, 
the value of L that we calculate will characterize the 
total system luminosity, not the luminosity of a single 
star. That is, unresolved binaries result in overestimates 
of the luminosity of the primary star. Since the value of 
sin i depends inversely on vL, the presence of unresolved 
binaries in our observational sample will cause sini to 
be systematically underestimated, or conversely, for the 
value of (sini) m to be systematically overestimated. 

To correct (sin?) m to account for unresolved binaries, 
we assume that the masses and luminosities of both the 
primary and the secondary can be related through a 
mass-luminosity relation of the form 



L = kM a , 



(7) 



where k and a are constants. Thus, we can express the 
total luminosity of the primary and secondary, Lt, as 



L T = L P (l + q a ), 



(8) 



where Lp is the luminosity of the primary and q is 
the mass ratio of the secondary to the primary, q — 
M s /M p < 1. Assuming that the temperature derived 
for the binary system corresponds to the temperature 
of the primary, then the value of (sinz) m calculated for 
the primary star in a binary system can be corrected by 
simply dividing by a correction factor of yT + q a . Since 
we do not know the exact values of q, we simply assume 
that q is uniformly distributed between and 1. While 
the true distribution of q is almost certainly not uniform, 
the deviation from a uniform distribution is likely small 
enough to not significantly impact our results. 

To determine the value of a in Eq. [8] we fit to the mass 
and luminosity models derived for pre- main sequence 
stars between the ages of 1-10 Myr by iBaraffe et al.1 
om this fit we derive a value of a = 1.5. 

The binary fraction for NGC 2264 is poorly con- 
strained. The multiplicity of nearby main-sequence 



field stars appears to be mass dependent, ran ging from 
~50% for G stars dDuauennoy fc Mavorlll991h to ~30% 
for M stars dFischer fc MarcvHl992t iReid fc GizidH997t 
iDelfosse et al.l l2004ft. The multiplicity for regions of iso- 
lated star formation, such as Taurus, has been found to 
be considerably highe r than the field p opulation, per- 
haps larger than 80% (|Simon et al.l ll993'l . Clustered re- 
gions such as the ONC and NGC 2264, however, show no 
such excess (|Kohler et al.l 120061 ). suggesting initial con- 
ditions or dynamical evolution at early ages have an im- 
portant effect of stellar multiplicity in the pre-main se- 
quence phase. As NGC 2264 is closer in character to the 
ONC than Taurus, we assume a value of B — 0.50 as the 
preferred binary fraction for our model, and investigate 
values as low as B = 0.0 and as high as B = 0.75. 

4.5. Comparing Model sini Distributions with 
Observations 

4.5.1. Measuring a best fit distance 

We compare the modeled distribution of sin i to the ob- 
served distribution using two-sided KS tests over a range 
of assumed distances. The KS tests were repeated 100 
times at distances of 600 pc to 1100 pc in steps of 5 pc. 
The distance with the highest median KS probability, 
Pks, is the best fit distance. Fig. [9] shows the Pks vs. 
distance curve for the most likely set of adopted model 
parameters (solid line: B = 0.5, all errors drawn from 
the observed error distributions). We find the best fit 
distance to be 913 pc. At this distance, the median KS 
probability is ~ 0.5. The observed and modeled distri- 
butions corresponding to this best fit distance are shown 
in Fig. 1 
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Fig. 8. — Comparison of the observed distribution of sin i assum- 
ing a distance to NGC 2264 of 913 pc with the modeled distribu- 
tion. 



4.5.2. Measuring statistical errors with a bootstrap analysis 

We have performed a bootstrap analysis to obtain a 
rigorous estimate of the statistical uncertainties associ- 
ated with our measurement of the distance to NGC 2264. 
In the bootstrap procedure, many artificial sini datasets 
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Fig. 9. — Effect of various assumed error models on the Pks 
vs. distance curve. Vertical axis represents the probability that 
the observed sini distribution at a particular distance was drawn 
from the same distribution as the modeled sin i distribution. 

are generated by selecting with replacement from the ac- 
tual set of observed sin i values. Each of these synthetic 
data sets are analyzed using the method described in 
i j4.5l to determine a best-fit distance to the cluster. The 
resultant distribution of derived distances provides an 
estimate of the true distance, and the width of the dis- 
tribution provides an estimate for the uncertainty in our 
derived distance. By using replacement, the bootstrap 
procedure tests the uncertainty in our distance estimate 
due to sampling effects, but does not account for any se- 
lection effects or biases that would introduce systematic 
differences between the true and observed sini distribu- 
tions. 

For the bootstrap analysis, we generated 200 sets of 
sini values of equal size to the observed sample (i.e. 
97 stars) by randomly selecting values with replacement 
from the list of observed sinz values. For each of the 
200 sets, a best fit distance was derived using the KS 
procedure described above. The resultant distribution of 
distance estimates is shown in Fig. [TO] We derive a l-cr 
confidence range for our distance estimate from the width 
of the velocity region that encloses 67% of the best-fit 
distances measured with the bootstrap procedure. Com- 
bining this uncertainty estimate with the best fit distance 
measured in £|4.5l vields a distance estimate to NGC 2264 
of 913 ±40 pc. 

4.5.3. Investigating systematic effects 

Our distance measurement includes statistical uncer- 
tainties due to sampling effects, as calculated above, but 
also potential systematic errors due to the assumptions 
that underlie our models. Factors that could introduce 
systematic errors into our analysis include: the error pre- 
scriptions and binary fraction adopted in our calculation 
of (sini) m ; the usini cutoff we imposed on our mod- 
eled and observed sini distributions; biases in the stellar 
properties derived for stars as a function of their evolu- 
tionary state; and finally, the underlying assumption of 
isotropically distributed rotation axes. We consider in 
turn the potential impact of each of these effects on our 
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Fig. 10. — The distribution of bootstrapped best-fit distances is 
overlaid on a plot of the Pks vs - distance curve for the observed 
data. A binary fraction of B = 0.5 is assumed, and errors are 
drawn randomly from the observed error distributions. 



analysis. 

The impact of different sources of observational error 
on the distance determination can be seen in Fig. [SJ 
where we compare the Pks vs. distance relations pro- 
duced by comparing our observed sample to models that 
neglect various components of the sin i error budget. The 
solid line represents the case where all of the errors are 
chosen randomly from the actual error distributions, and 
the remaining curves represent cases where the errors in 
a particular variable have been set to 0. It is clear from 
Fig. [9] that the assumed errors have a non- negligible im- 
pact on both the best-fit distance as well as the height 
of the Pks curve. For instance, negating the luminos- 
ity errors has the effect of reducing the best-fit distance 
by 4 %, while eliminating the T e // errors increases the 
best-fit distance by roughly G %. Eliminating the usim 
errors not only increases the best fit distance by 4 %, but 
also strongly reduces the peak Pks, suggesting that the 
errors in vs'mi contribute significantly to the shape of 
the sin i distribution. 

To the extent that the error distributions we adopt in 
our model (see ^4. 3|) do not reflect the true error distri- 
butions of our data, our distance estimate will be skewed. 
As shown above, the error distributions we do adopt are 
influencing our derived distance at the ~5% level, sug- 
gesting any systematic error in our derived distance due 
to adopting improper error distributions would likely be 
~5% as well. Our error distributions, however, are con- 
sistent with those estimated by other authors for the 
same parameters, and the error distributions push the 
derived distance in different directions, such that multi- 
ple systematic errors should offset one another to some 
degree. Perhaps most worrisome is the asymmetric shape 
of the luminosity error distribution, but lacking a more 
robust means of characterizing the luminosity errors, we 
are unable to remove this potential systematic effect from 
our analysis. 

Fig. HU shows the effect of the assumed binary frac- 
tion on the best fit distance. As the binary fraction is 



increased, the average luminosity of the modeled stars is 
increased, causing the sinz distribution to move to lower 
sin i as per Eq. [5] In order for the observed sin i distribu- 
tion to remain well matched to the modeled distribution, 
the assumed distance must therefore increase according 
to Eq. [5] Thus, we see that as the modeled binary 
fraction changes from B — 0.0 to B = 0.75, the best- 
fit distance goes from ~ 860 pc to ~ 970 pc, a change 
of roug h ly 13 %. This is consistent with the results of 
IJeffriesI (|2007| ). who found in his ONC study that in- 
creasing the binary fraction by 0.2 led to a 4% increase 
in his modeled best fit distance. 
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Fig. 11. — Effect of various assumed binary fractions, B, on the 
PkS vs - distance curve. Vertical axis represents the probability 
that the observed sini distribution at a particular distance was 
drawn from the same distribution as the modeled sini distribution. 



Biases in the parameters derived for stars in different 
evolutionary stat es could also t he results of our analy- 
sis. As noted by IJeffriesI (|2007), it may be more diffi- 
cult to derive accurate temperatures, luminosities, and, 
thus, radii for classical T Tauri stars than more evolved 
weak T Tauri stars; classical T Tauri stars typically pos- 
sess larger extinctions, are irregularly variable, and have 
significant contributions to their total luminosity from 
accretion and re-radiation from their circumstellar disk, 
all of which can complicate the der ivation of their in - 
trinsic stellar properties. Following lLada et al.l (2006), 
we used an IRAC [3.6 - 8] color cut to identify classi- 
cal T Tauri stars in our sample; when these CTTSs are 
excluded from our sample, the best fit distance to NGC 
2264 increas es to ~950 pc . This effect is contrary to 
that seen by I Jeffrie"sl (|2007l ). who found a decreased best 
fit distance to the ONC once CTTSs were excluded from 
his sample. Given these contradictory results, and that 
the effect of removing CTTSs from our sample is compa- 
rable to the other systematic effects probed here, we have 
chosen not to exclude CTTSs from our main analysis. 

Our technique assumes that all of the stars in NGC 
2264 are at the same distance from Earth. In reality, the 
cluster has some li ne of sight d epth, predicted to be on 
the order of 28 pc (|Dahml 12008). If we assume that the 
known members of NGC 2264 are biased towards low ex- 
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tinctions, this could cause the distance we derive to the 
cluster to be slightly smaller than the true geometric cen- 
ter of the cluster. The distance we derive here, however, 
would properly describe the distance of the population of 
currently known members; we do not consider this to be 
a systematic error in our analysis, but rather a nuance 
that should inform the interpretation of our results. 

Finally, the fundamental assumption of random axial 
orientations is itself somewhat questionable since stellar 
clusters like NGC 2264 are predicted to have collapsed 
from single cloud complexes. We might expect, then, 
for there to be a preferred orientation of stars in the 
cluster resulting from the conservation of the cloud's ini- 
tial angular momentum, or possibly from the presence 
of large scale magnetic fields. Previous applications of 
the sini distribution technique have, however, produced 
results that are in agreement with precise paralla x mea- 
surements (|Preibisch fc Smithlll997l Ijeffriesll2007h . This 
agreement provides evidence in support of the random 
axial orientation assumption, but it has not been directly 
confirmed by observations. 

The discussion above has revealed a number of poten- 
tial systematic effects in our analysis. The five potential 
systematics most amenable to direct investigation (the 
inclusion of CTTSs in our primary sample; the adopted 
luminosity, temperature, and usini error distributions; 
and the assumed binary fraction) all influence the de- 
rived best fit distance at the ~ 5% level. We therefore 
combine these individual uncertanties in quadrature to 
characterize the potential error in our best fit distance 
due to systematic effects, producing a estimate of our 
total systematic uncertainties of ± 12 %, or 110 pc. 

4.6. The derived distance and age of NGC 2264 

We have calculated distance-dependent sin i values for 
a sample of 97 pre-main sequence stars in the open clus- 
ter NGC 2264. By comparing the observed sini distri- 
bution to a modeled distribution assuming random ax- 
ial orientations, we derive a distance of 913 pc to NGC 
2264; quantitative tests of our analysis reveal sampling 
and systematic errors of 40 and 110 pc, respectively. Our 
distance estimate does not rely on evolutionary models 
to any significant degree. 4 

Our distance estimate is significantly higher than a 
number of previously determined distances, particularly 
the w idely cited value of 760 pc found bv ISung et al.1 
(|1997f ). In general, though, our estimate falls within the 
typical range of calculated distances (730 pc to 950 pc) 
for NGC 2264. Table [2] provides a comparison of our re- 
sult to previous distance estimates. A distance to NGC 
2264 of 913 pc represents an increase of approximately 
20% compared to the widely excepted value of 760 pc, 
though the two results are formally consistent within the 
sum of the statistical and systematic error bars. 

The mean age of NGC 2264 is commonly cited as ~ 
3 Myrs, though there is ev idence for a considerable age 
spread within the cluster (|Dahml [2008D . The luminosi- 

4 The one exception is the use of the pre-main sequence mass- 
luminosity relation in the derivation of the binary correction factor. 
This factor had a relatively small effect on our calculated distance, 



ties of pre-main sequence stars are often compared with 
predictions of theoretical pre-main sequence models to 
infer the age of their parent cluster; as luminosity de- 
clines through the pre-main sequence phase, the larger 
luminosities produced by assuming a greater distance to 
the cluster will produce a younger inferred age for the 
cluster. Increasing the assumed distance to NGC 2264 
from 760 to 910 pc changes the distance modulus by 0.4 
mag. The corresponding 0.4 mag brightening of the stars 
produces a shift in the age of the cluster. We have pro- 
duced a crude estimate of the revised age of the clus- 
ter by determining the age at which a 1 M© star's H 
band magnitude is 0.4 mag brighter than at 3 Myrs: ac- 
cording to the pre- main sequence models calculated by 
iBaraffe et al.l (|1998t ) , the distance we derive here implies 
a revised age for NGC 2264 of -1.5 Myrs. The detailed 
analysis of cluster members required for robust estimates 
of the age and properties of NGC 2264 in light of our 
new derived distance, however, is beyond the scope of 
this work. 

5. CONCLUSIONS 

We determined the distance to the open cluster NGC 
2264 using a statistical analysis of cluster member in- 
clinations. We derived distance-dependent values of sini 
(where i is the inclination angle) for 97 stars in NGC 2264 
from measured rotation periods, luminosities, effective 
temperatures, and projected equatorial rotation veloci- 
ties, i»sini, of these stars. We measured 96 of the vsini 
values in our sample by analyzing high-resolution spec- 
tra with a cross-correlation technique. We modeled the 
observed distribution of sin i for the cluster by assuming 
that member stars have random axial orientations and 
by adopting prescriptions for the measurement errors in 
our sample. By adjusting the distance assumed in the 
observed sini distribution until it matches the modeled 
distribution, we obtained a best-fit distance for the clus- 
ter. We find the data to be consistent with a distance to 
NGC 2264 of 913 pc. Quantitative tests of our analysis 
reveals uncertainties of 40 and 110 pc due to sampling 
and systematic effects, respectively. This distance esti- 
mate suggests a revised age for the cluster of —1.5 Myrs, 
although more detailed investigations of the full cluster 
membership are required to draw strong conclusions. 

The authors wish to thank Steve Strom for a prompt 
and helpful referee report that improved the analysis pre- 
sented here, and Russel White and Jeff Burchfield for 
providing IDL code that formed the basis of our cross 
correlation pipeline. EJB acknowledges the support of 
the SAO Summer Intern Program, made possible by a 
grant from the NSF. NASA support was provided to K. 
Covey for this work through the Spitzer Space Telescope 
Fellowship Program, through a contract issued by the 
Jet Propulsion Laboratory, California Institute of Tech- 
nology under a contract with NASA. 

however, and is also well constrained by observations. Thus, our 
derived distance is almost entirely model-independent. 
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TABLE 3 

Data used to create sini distribution 



RA 

(dcg) 


DEC 

(deg) 


log L bo i 
(ergs/sec) 


\ogT eff 
(K) 


Lboi and T e // 
Source 


Period 
(days) 


v sin i 
(km/s) 


tjsini Error 
(km/s) 


99.8594131 


9.6863384 


33.590 


3.630 


S 


5.51 


16.5 


1.51 


99.8765793 


9.5604086 


33.530 


3.640 


s 


5.49 


57.5 


6.39 


99.9138336 


9.9332304 


33.420 


3.590 


S 


0.86 


66.9 


5.33 


99.9232178 


9.5779305 


34.180 
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